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Henkel and Pleimling reply: In the preceeding 
comment, Corberi, Lippiello and Zannetti [1] studied 
the scahng of the hnear autoresponse function s) = 
s~'^~°'fRit/s) of ageing systems in the scahng regime of a 
quench below Tc- Here s is the waiting time, t the obser- 
vation time, fa a scaling function and a an exponent to 
be determined. The exponent a also describes the scaling 
of the thermoremanent magnetization 


TRM 


/h 


duR{t,u) - s-Tit/s) 


(1) 


where h is the external magnetic field. While we had 
found previously that scaling is observed if a = 1/2 is 
taken [2], they assert that the asymptotic value of a 
should be different. Specifically, they suggest a = 1/4 
and a = 1/2 with logarithmic corrections, for the 2D 
and 3D kinetic Ising model, respectively [1]. This asser- 
tion is supported by numerical data for the field-cooled 
magnetization MzFc/h = J^du R{t,u) ^ s^°'G{t/s) ob- 
tained for both the 2D and 3D Ising model. In addition, 
they argue that for Mtrm, the true asymptotic scaling 
behaviour should set in only for very large values of s. 
Instead, they suggest that their data rather fall into a 
preasymptotic regime with an effective Oeff = A/z, where 
z and A are the dynamical and Fisher-Huse exponents [3] , 
respectively. This is supported by numerical data going 
up to waiting times s ~ 200 in 3D and s ~ 2 • 10'^ in 2D. 
The truly asymptotic regime is not yet reached by their 
data, but they expect to recover the same value of a as 
obtained from Mzfc [!]■ 

In order to test the proposition of [1], we obtained data 
for Mtrm for the 2D Ising model at zero temperature, for 
very long waiting times up to s = 5600. In analysing the 
scaling of Mtrm, we did not perform any subtraction 
in order to isolate an 'ageing part', in contrast to the 
procedure adopted in [1]. Our result is shown in Figure 1, 
for several values of the scaling variable x = t/s. To 
guide the eye, we also show the power laws s"^/'', s""**/^ 
(A/z = 0.625 in 2D [3]) and fi-^/^. In figure 2 we focus 
on the scaling behaviour for large values of s. 

At first sight, the data from Figure 1 appear to be 
roughly consistent with a simple power law (1), with 
a between 0.5 and 0.625. On closer inspection, three 
regimes are visible. In the first regime, for small waiting 
times, an approximate scaling with an effective Oog ~ 1/2 
is obtained. A second well-defined scaling regime ex- 
ists for intermediate times (here approximately between 


10 < s < lO'^) and Mtrm ~ s" 


Finally, for very 


long times, there is a cross-over into a third regime, where 
TRM ^ see figure 2. The same scaling regimes 


M' 


can be found for any temperature T < Tc [5], as expected, 
see [4] , but data for T > are much more difficult to ob- 
tain. We do not find any hint for a regime with a = 1/4. 
In [1], working with smaller values of s, only the first two 
regimes were observed. 


We also point out that the functional form of Jr^x) 
(and consequently F{x)) is independent of s. 

In conclusion, we find evidence that the truly asymp- 
totic scaling of the thermoremanent magnetization 
should be described by an exponent a — 1/2 in the 2D 
Ising model. This is different from the value suggested 
from the analysis of Mzfc [!]■ This work was supported 
by CINES Montpellier (projet pmn2095). 
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FIG. 1. Scaling of Mtrm as a function of s for the 
zero-temperature 2D Ising model. 
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FIG. 2. Power-law scaling of Mtrm for large values of s. 
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